Phase dynamics in a binary-collisions atom laser scheme 



0. Zobay and P. Meystre 
Optical Sciences Center, University of Arizona, Tucson, AZ 85721 



OS 

o> 

> 
o 

(N 
(N 



m 

(N 



o> 

c3 



c 

o 
o 



Various aspects of the phase dynamics of an atom laser 
scheme based on binary collisions are investigated. Analytical 
estimates of the influence of elastic atom-atom collisions on 
the laser linewidth are given, and linewidths achievable in a 
recently proposed atom laser scheme [Phys. Rev. A 56, 2989 
(1997)] are evaluated explicitly. The extent to which a relative 
phase can be established between two interfering atom lasers, 
as well as the properties of that phase, are also investigated. 

PACS numbers: 42.50.Vk,03.75.Fi,32.80.Pj,42.50.Ct 



I. INTRODUCTION 

The recent advances in the generation and manipula- 
tion of samples of ultracold atoms have stimulated re- 
search aimed at the development of sources of coherent 



atomic beams. Recent MIT experiments [ljQ] can be con- 
sidered as the first demonstration of a pulsed atom laser. 
An important further goal is to achieve a cw system. Sev- 
eral theoretical proposals to realize such a device have 
already been made. They can be separated into two cat- 
egories: in the first one, the population of the atom laser 
mode results from optical cooling-like processes such that 
electronically excited atoms undergo a transition to the 
atom laser mode via spontaneous emission || -|?J . Pro- 
posals in the second category |§Hl2]] rely on binary colli- 
sions between ground state atoms fed into an intermedi- 
ary level of the atom resonator. These collisions are such 
that one of the atoms undergoes a transition into the laser 
mode and the other into a heavily damped mode. This 
latter atom rapidly escapes the cavity, thereby providing 
the irreversibility of the pumping process. 

Amongst of the most important characteristics of a 
laser are the coherence properties of its output. In the 
context of cw atom lasers, this aspect has been addressed 
in a number of articles: Refs. ]l^-|lf|| investigate vari- 
ous output coupling schemes and their effect on the laser 
linewidth. They show that a narrow linewidth can be 
achieved by a suitable choice of outcoupling setup. How- 
ever, the details of the internal dynamics of the atom 
laser have not been taken into account so far in that work. 
In contrast, Refs. [||jl(| consider the correlation functions 
for a three-mode model of a binary-collision atom laser, 
modeling the output mechanism simply as linear damp- 
ing. These studies show that the elastic atomic collisions 
which occur in addition to the pumping collisions impose 
severe limitations on the coherence of the laser output. 
Finally, estimates of the linewidth of an atom laser based 
on optical cooling are given in Ref. |]. 



The purpose of the present article is to extend the 
previous work on the phase dynamics of a binary col- 
lision atom laser. We describe this system in terms of 
a three-mode scheme which ignores most of the multi- 
tude of atomic cavity levels and takes into account only 
those modes that are essential for the lasing process, i. e. 
the pump, laser mode and loss mode. The feeding and 
loss processes are described by linear coupling to external 
continua. 

Several reasons motivate our interest in the binary col- 
lision atom laser and in this particular model. First, and 
in contrast to the optical pumping approach, evaporative 
cooling has already demonstrated its value in achieving 
BEC 17 19j . The binary collision model we consider can 
be seen as the simplest possible caricature of this mecha- 
nism. Second, an interesting proposal has recently been 
made to overcome some of the problems associated with 
elastic binary collisions |ll],[l2| in case they result from 
the dipolc-dipole interaction between atoms: By taking 
advantage of specifically designed optical cavities, it is 
possible to greatly decrease the rate of elastic collisions 
between laser mode atoms. It is of interest to examine 
the resulting influence on the laser linewidth. Third, the 
binary collision atom lasers can be analyzed within the 
framework of few-mode models in more detail and more 
easily than the models based on optical cooling. An addi- 
tional justification for studying such rudimentary models 
lies in the fact that they give a good theoretical under- 
standing of the fundamental way an atom laser might 
generically work. As such, they provide valuable com- 
plementary insights to more realistic approaches, e.g. the 
quantum kinetic theory of BEC which may also be ap- 



plied to atom lasers [20l 22 . Finally, before including the 



effects of more sophisticated output coupling schemes it is 
important to gain a thorough understanding of the inter- 
nal laser dynamics and its impact on the laser linewidth 
in simple models. 

After briefly reviewing the main aspects of binary col- 
lisions atom laser schemes in section II, the present paper 
addresses two main topics: Section III discusses the laser 
linewidth and its dependence on elastic atom-atom colli- 
sions. In previous studies this question was either treated 
rather briefly jl0| or estimates were given || the valid- 
ity of which is not quite clear [^3| . Our analysis leads to 
an analytical approximation for the laser linewidth in the 
presence of interatomic collisions. To this end a linearized 
fluctuation analysis of the system master equation is per- 
formed after adiabatically eliminating the loss mode. The 
results of this analysis are compared to quantum Monte 
Carlo simulations [ p4p^ ]. We illustrate these results in 
the determination of the possible operating regimes of 
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the atom laser of Refs. 0,Q, and demonstrate that in 
for weak pumping it should in principle be possible to ob- 
tain laser linewidths below the natural linewidth of the 
laser mode. 

In section IV we turn to the question of whether a 
definite relative phase can be established between inde- 
pendent atom lasers. The same question, but for Bose- 
Einstein condensates, has aroused much interest recently 
pG| pl{ . It was shown that for a system of two conden- 
sates with a fixed total number of atoms, a relative phase 
is established by the process of measurement. Our anal- 
ysis extends this work to the realm of open systems. We 
show that a relative phase can be established in single 
runs of the experiment in very much the same way as for 
Bose condensates. The study of the diffusion properties 
of this relative phase shows its close connection to the 
phase of the single atom laser as described in terms of 
stochastic processes. Finally, a brief summary and con- 
clusions are given in section V. 

II. THE BINARY COLLISION ATOM LASER 
MODEL 

The basic principles of operation of the binary collision 
atom laser are discussed in detail in several publications 
|^| fL2| |, One considers a resonator for atoms, realized e.g. 
by optical fields. In order to concentrate on the essen- 
tial dynamics only three out of the multitude of atomic 
center-of-mass modes are taken into account explicitly 
(cf. Fig. 0). Bosonic atoms in their ground electronic 
state are pumped into an atomic resonator level of "in- 
termediary" energy (mode 1). They then undergo bi- 
nary collisions which take one of the atoms involved to 
the tightly bound laser mode 0, whereas the other one 
is transferred to the heavily damped loss mode 2. This 
latter atom leaves the resonator quickly, thereby provid- 
ing the irreversibility of the pumping process. A macro- 
scopic population of the laser mode can built up as soon 
as the influx of atoms due to pumping compensates for 
the losses induced by the damping. 

In the description of this laser scheme one has to 
take into account that in addition to the pumping colli- 
sions other types of interatomic collisions can also occur. 
These considerations lead to an ansatz for the atom laser 
master equation of the form 

W = -i[H + H coi , W] + K V[a ]W + Kl (N+ l)V[a{\W 
+K 1 NV[a\]W + K 2 V[a 2 ]W (1) 

with h = 1. In this equation, we use the second quantized 
formalism in which each center-of-mass atomic mode is 
associated with an annihilation operator ai, and W de- 
notes the atomic density operatorn. The free Hamilto- 



1 Since we consider ground state atoms only, they are fully 



nian is given by 

H = ^2 Uia\ai, 

4=0,1,2 

u>i being the mode frequencies. The general form of the 
collision Hamiltonian can be written as 

Hcoi = V i:j kia\a]akai (2) 

i<j,k<l 

with Vijki the matrix elements of the two-body interac- 
tion Hamiltonian responsible for the collisions. However, 
for the present purposes it will be sufficient to restrict 
our attention to the reduced form 

H co i = V 2iialala 1 ai + V 110 2a\a{a Q a2 + Vboooao a o a o a o 
+Vbioia a i a o a i + Vima\a\aiai (3) 

in which (besides the pumping collisions) only those col- 
lisions are retained which are expected to have the most 
significant influence on the phase dynamics. The damp- 
ing rates of the cavity modes are given by the coefficients 
Ki, and the strength of the external pumping of mode I 
is characterized by the parameter N, which is the mean 
number of atoms to which mode 1 would equilibrate in 
the absence of collisions. The superoperator T> is of the 
Lindblad form and is defined by 

V[c]P = cPc) - \ (c+cP - Pet c) (4) 

with arbitrary operators c and P. 

It has been pointed out || that it is useful to study this 
system in detail although it might be thought of more as 
a caricature than as an approximation, due to the neglect 
of most of the modes in the atomic resonator. From this 
point of view, its main purpose is to give a basic un- 
derstanding of the laser dynamics, rather than a detailed 
quantitative description of a realistic setup. However, re- 
cent work has added significantly more substance to this 
model [ pd|JT^ ]: It was shown that in case the binary col- 
lisions result from the dipole-dipole interaction and for 
a quasi one-dimensional modulated cavity, a high mode 
selectivity for the atomic collisions can be obtained. One 
particular pumping mode is then coupled only to a small 
number of other modes. Among the small set of initially 
populated low-lying modes, a particular one is then ex- 
pected to be singled out as a result of mode competition. 
Another benefit of this model is the fact that due to the 
dipole-dipole selection rules in the resonator, the strength 
Voooo of detrimental collisions can be made small in com- 
parison to Vo2ii- In that proposal, an efficient pump- 
ing process is achieved by a time-dependent modulation 
of the atom cavity. In the spirit of this model we set 
<-^o + ^2 = 2t^i in the following. 



described by their center-of-mass quantum numbers. 
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In order to achieve a sufficiently high degree of irre- 
versibility it is necessary that Ki is much larger than the 
damping rates of the other modes. This suggests to adi- 
abatically eliminate this mode, an approximation that 
leads to the simplified master equation |g|-|Tc[] 



-i[H c ,p] + K V[ao]p + Ki(N ■ 
-KiNV[a\]p + TV[ala\]p. 



l)V[ ai ]p 



(5) 



Equation (0) is written in the interaction picture with 
respect to Hq = wofigao+^i&ifli; and the reduced density 
matrix p is p — Tr moc j e 2^]- The reduced collision 
Hamiltonian H c is 

H c = Vbooo<4 a o a o a o + Vbioiao a i a o a i + Vimajajaiai, 

(6) 

and T = 4\V 2ii\ 2 /k. Consistently with Ref. § we call 
the limiting cases r <C kq and r 3> kq the weak and 
strong collision regimes, respectively. The reduced mas- 
ter equation (^|) forms the basis of most previous studies 
of binary collision atom lasers. 

The master equations ([!]) and (g) can easily be solved 
numerically using standard quantum Monte Carlo sim- 
ulation techniques Q,^5|. This is discussed in detail in 
Ref. 1 12 1, and we only recall that these simulations are 
significantly facilitated by the fact that between quan- 
tum jumps the total number of atoms is conserved. In 
this way one has to propagate only a small number of 
actually populated states, a situation evidently advan- 
tageous with regard to both memory and CPU time re- 
quirements^]. 



III. ELASTIC COLLISIONS AND LASER 
LINEWIDTH 

A. Linearized fluctuation analysis for the two-mode 

system 

In order to obtain an analytical approximation for 
the laser linewidth in the two-mode sytem a linearized 



fluctuation analysis can be performed f p2|j33| ]. To this 
end the master equation (||) is converted to a Fokker- 
Planck equation using the P-function representation as 
described in |1. This equation can be transformed to 
polar coordinates ctj — y/rTje l ^ j , where otj denotes the 
complex amplitudes originally appearing in the Fokker- 
Planck equation |34|,|35|]. This leads to stochastic differ- 
ential equations 

dn = [Tnl(n a + 1) - n n Q ]dt + dS no , (7) 

dm = [Ki(N-m) - 2Tnl{n + l)]dt + dS nil (8) 

#o = [-Vbooo(2no - 1) ~ Voioini]dt + dS</> , (9) 

c% = [-V r im(2ni-l)-Vbioino]* + ^ 1 . (10) 



The correlation matrix for the stochastic forces dS T = 
'-Aoj^Vi) is given by 



(dS no , dS ni , dS& 



D 



/ 2Tn{n -2Tnln 
-2Tn\n 2K 1 Nni - 2Tn\n 
-2Vbooo"o -VqwiUi 
\ -Voioi^o — 2V\u\ni 



-2V ooo"-o -Vbioi^o \ 
-Vbioi^i -2T4mni 
Tnl/(2n ) Trn/2 

rni/2 ^ + V -T j 
(11) 



In the limit no ^ 1 one obtains from Eqs. (Q) and 
(||) the above-threshold semiclassical steady-state pop- 
ulations |o) 



n 



ni 




(12) 
(13) 



the threshold condition being TV > ^Jkq/T. The re- 
lation (|lj) between no and fi\ also holds in the full 
quantum- mechanical two- and three-mode models [ p"2| . 
To proceed further we introduce the fluctuation vari- 



ables Si 



and 



where the 



phase drift variables <f>j obey the deterministic equa- 
tions obtained from Eqs. (||) and ( |l0| ) by discarding the 
stochastic forces and substituting hj for rij. In the lin- 
ear approximation the time evolution of the fluctuations 
SA T = (5no, 5ni, 54>o, S4>i) is given by 



dSA = -kSAdt + dS 



(14) 



In particular, in the calculation of correlation functions of 
the type (aUt + T)aj(t)), one normally has to evolve four wave 



where the matrix k is obtained by linearizing the drift 
terms in Eqs. (0) - ([l0]) around the steady-state values 
The correlation matrix for the stochastic forces dS is 



functions [l±(i) 0j ]|$(t)) for a given initial sample wave func- giyen by the matrix D of Eq (|fj|) after replacing rij by 



tion |3>(t)} if one uses the method of Ref. J2f|]. However, due 
to the conservation of atom number in this case the contri- 
butions from |4>(t)) and a 3 |$(i)) in the time development of 
[1 + a,-]|<&(t)) can be distinguished. This means that one ob- 
tains a contribution to the correlation function even from a 
single propagation. As it is more advantageous to average over 
a large number of different |3?(t)) than to use many simula- 
tions for a small number of initial wave functions one obtains 
a significant additional decrease in computation time. 



Supposing that SA — at time r = 0, the probability 
distribution of the fluctuations at a later time r is given 

by a 



p(8A,t) = [(27r) 4 det£r(r)]- 1/2 exp[ 



with 



■\5A T a- 1 { T )5A] 
(15) 
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ct(t) = / dr' exp[-k(r - r')] D cxp[-k T (r - r')]. 
Jo 

(16) 

Assuming as usual that the first-order correlation func- 
tion Cq(t) = (aJ(r)ao(O)) is determined only by the 
phase fluctuations one obtains 



C (t) = I d 4 5An e-^°p(5A,r) 
= n e"* exp[-|«733(T)], 



(17) 



where the deterministic phase drift is also included. This 
phase drift leads to a shift of the center of the power spec- 
trum (the Fourier transform of the correlation function) 
by an amount 2^000^0 + Voioi^i with respect to the col- 
lisionless case. The behavior of the correlation function 
Cq(t) is thus essentially determined by the covariancc 
matrix element 1733 (r). 

The right hand side of Eqs. (0) - ( |i~Ci| ) depends only on 
the atom numbers no and niR This means that in that 
two-mode scheme the phase diffusion due to interatomic 
collisions is induced solely by atom number fluctuations. 
This observation significantly simplifies the further ana- 
lytical treatment as it leads to the matrix k having two 
vanishing column vectors. Hence the covariance matrix 
<t(t) can be evaluated explicitly in a straightforward way 
according to Eq. However, the ensuing expression 

is still rather complicated, so that in the following we 
restrict the discussion to two limiting cases which illus- 
trate the essential aspects of the influence of the elastic 
collisions on the laser linewidth. It should also be noted 
at this point that in the two-mode system Co(r) does 
not depend on elastic collisions between pumping mode 
atoms, which are characterized by the parameter Vim. 

(a) V0101 = 0. Expanding the expression for 033(7") to 
leading order in the parameter no one obtains 



0-33(7") = t[w + k / (2n )] 
[1 - cxp(-qr)] 



(8V 2 000 n /g- 2w) 



with 



w 



and 



V 2 



[1 - exp(-2gr)] 
2q 



KlN 



(w - 8V 2 000 n /q) (18) 



2 + 2yWf- 



N 



(19) 



3 This is in contrast to the three- mode system, where an 
explicit phase dependence is introduced through the coherent 
pumping term Vo2iiOo a 2 a i a i + ^- c - 



Q : 



44 



4k + Ki"-i/n 



(20) 



^From Eq. ( |l2| ) it follows that KiN and thus w are of 
the order of uq. In the first line of Eq. ( |l8| ) the term 
Ko/(2n,o)r was included although it is not of the same 
order in no as the other terms. This is because it de- 
scribes the behavior of the correlation function if V0000 
is very small (cf. Fig. [|). The parameter q of Eq. ( |20| ) 
is one of the eigenvalues of the upper left 2 x 2-minor of 
the matrix k. It can thus be interpreted as the inverse 
of one of the timescales relevant in the dynamics of the 
atom number fluctuations. The other eigenvalue is of the 
order of no and introduces a much shorter time scale. It 
does not play a significant role for the characterization 
of the correlation function. 

Equation ( |l7j ) implies that the most important aspects 
of the correlation function can be inferred from the study 
of the behavior of 033(7") in the time interval where it 
is smaller than or of the order of unity. From Eq. dl§| ) 
two different kinds of behavior can thus be distinguished, 
depending on the magnitude of Voooo ■ As long as 033 (r = 
1 /q) <C 1 the time evolution of (733 (t) relevant for Co (t) 
is well approximated by 



0-33(7") ~ t[w + K /(2n )]. 



(21) 



This is because as time increases the second and the third 
term on the right-hand side of Eq. ( |ilf ) remain constant 
for r > 1/q, while the first term increases. For such 
values of Voooo Co (r) decays therefore exponentially, and 
the power spectrum is Lorentzian. If w S> fto/(2no) the 
linewidth is proportional to V ooo"-o- IN contrast to the 
situation with conventional lasers, it increases linearly 
with the number of atoms in the laser mode. 

In case <t 33 (t = 1/q) 3> 1 the decay of the correlation 
function can accurately be approximated by expanding 
the exponentials in Eq. (|l^) up to second order. This 
yields the expression 



0-33(7") ^ 4V 2 000 n T 2 



(22) 



Under these circumstances the correlation function de- 
cays like a Gaussian. The spectrum is thus itself a Gaus- 
sian and its linewidth is proportional to Voooo y/no- The 
atom laser linewidth still increases with no, albeit less 
dramatically than in the preceding case. 

Figure || suggests that a rough estimate of the value 
of Voooo at which the decay of the correlation function 
changes from exponential to Gaussian can be obtained 
in the following way. One determines Voooo such that 
both Eqs. (21) and ( p2| ) yield the same solution for the 
condition 0-33 (r) = 2/ In 2. In this way one obtains as 
"critical" value 



V 



0000, crit 



/n 



(23) 



with w' = w/ V ooo- ^ n ° 1S mcrea sed the change thus 
occurs for smaller values of V 



0000- 
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It is instructive to compare these results to those ob- 
tained for a damped harmonic oscillator with Hamilto- 
nian H = ujoo^clq + Voooo«o a o a o a o- Denoting in analogy 
with Eq. ([!]) the damping coefficient by kq and the ex- 
ternal pumping strength by N , a linearized fluctuations 
analysis yields for the relevant covariance matrix element 



a(r) 



RV 2 N* + N 

O Vr 



'0000 



K 



.16V 3 - 2+N/2 

iov 0000 



N 



k /2 

-, [1 - exp(-K r)] 



+ 4 K5boo-r[ 1 - cx P(- 2K or)] 



(24) 



o 



where no further approximation has been made, in con- 
trast to the case of Eq. (|l8|). In the limit of large Vbooo 
Eq. yields ct(t) = AV^ Nt 2 , a result formally iden- 
tical to Eq. ( |22] ) since the equilibrium population uq be- 
comes N in the one-mode system. Such a connection 
between the linewidths for this one-mode system and the 
atom laser has been previously noted in Ref . p3j . Despite 
their formal similarity, there are however significant dif- 
ferences between Eqs. (|l8| ) and (|I|). First, there is no ob- 
vious correspondence between the coefficients multiply- 
ing their respective time-dependent contributions, and 
they are not of the same order in uq. In addition, the 
one-mode system is characterized by a very large "funda- 
mental" linewidth kq/2 which masks a possible quadratic 
dependence on Voooo- Finally, the above-mentioned lim- 
iting behavior is reached for values of Vbooo much larger 
than in the two-mode atom laser system. 

(b) Vbooo = 0. In this case the expansion of 0-33 (r) to 
leading order in uq yields 



^33(7) 



( ^0101 , K o 
\ 2Tno 2fiQ 



(25) 



The correlation function thus decays exponentially for 
all values of Voioi- A qualitative change in behavior as 
in the previous situation does not occur. The linewidth 
of the spectrum is now proportional to Vo 2 10 i/^o, an d 
becomes narrower when the population of the laser mode 
is increased, very much like the familiar Shawlow-Towncs 
linewidth of conventional lasers. 

The analytical results described so far were compared 
to numerical calculations of the correlation function. 
Both the weak-collision regime and the strong-collision 
regime were investigated for the explicit values T/kq = 
1/15 and 15 and K\/kq = 10 and 100. For Vbioi = a 
good agreement could generally be observed between the 
numerical result and the analytical approximation based 
on Eq. ( ]l7| ) . A representative example is shown in Fig. ^ 
where the parameter values are T/ko = 15, ki/kq = 100, 
N = 1.85, Voooo/ K o = 2.5, V0101 = 0. The numerical 
results are well approximated by a Gaussian. Consider- 
able deviations occurred only in the case T/ko = 1/15, 
ki/kq = 100. For these values the quantum effects in the 



population dynamics are particularly large, e.g. choosing 
N such that fio = 50 according to Eq. ( |l2| ) led to a nu- 
merical value of 100 whereas in the other case deviations 
were only a few percent. It should be noted that in the 
two-mode systen the population dynamics do not depend 
on the interatomic collisions. 

The transition between the exponential and Gaussian 
regimes in the behavior of the correlation function Co (r) 
is illustrated in Fig. |[ It shows the inverse of the "half- 
life time" T X [ 2 of Co(r) (defined by 0-33(^/2) = In 2/2) 
as a function of Voooo for the parameters T/kq = 15, 
K i/ K o = 10, and N = 20.3, corresponding to no = 100, 
on a doubly-logarithmic scale. The full curve depicts Ty 2 
as determined from Eq. (|l7|) , the dashed and the dotted 
curves are obtained from the approximate relations j2l| ) 
and (p2]), respectively. The results of numerical Monte 
Carlo simulations are shown as circles (•). The figure in- 
dicates that the analytical predictions, in particular the 
existence of two different regimes and the transition be- 
tween them, are accurately confirmed by the numerical 
calculations. The approximations based on Eqs. (^lj) and 
(P2I ) are almost indistinguishable from the full analytical 
expression in their respective regions of validity. 

The results of a similar study for the collision coeffi- 
cient Vbioi are shown in Fig. f|. There, the dashed curve is 
derived from the approximation (25). The global behav- 
ior of the numerical results is well approximated by the 
analytical description, in particular the quadratic depen- 
dence of the half-life time on Vbioi is recovered. However, 
the analytical prediction typically overestimates the nu- 
merical half-life time by a factor of two. This behavior 
is also observed in examples with other parameter val- 
ues. This indicates that the atom number fluctuations in 
the pump mode are not as well described by a linearized 
ansatz as the fluctuations in the laser mode. 

In case both collision coefficients Voooo an d V0101 are 
non-vanishing the analytical form of a reliable approxi- 
mation to Eq. ([l6]) would be rather complicated. How- 
ever, it can be seen that all terms involving Vbioi are at 
most of the order of Uq. Taking into account Eq. (^3|) this 
means that in the limit of large hq the function 033(1") 
can eventually be approximated by Eq. (p2|). The com- 
parison with numerical results shows that Eq. ( |l7| ) is still 
accurate when both collision coefficients are present. The 
calculations indicate that the correlation time decreases, 
i.e., the linewidth broadens, if the second collision coeffi- 
cient is switched on in addition to the first. If the values 
of the collision coefficients are such that one of them has 
a much larger impact on the correlation time than the 
other - which is checked by comparing Eqs. ( |l§| ) and 
( p5| ) - then the approximate results of the limiting cases 
(a) and (b) still hold. If both have roughly the same im- 
pact then the ensuing correlation time will still be of the 
order of the result inferred from the approximations. 
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B. Operating regimes of the atom laser 

In the following we apply these results to the atom 
laser scheme of Refs. [p p|Jl^ ]. In this model the col- 
lision mechanism is the laser-induced dipole-dipole in- 
teraction. Due to energy and momentum conservation, 
the associated selection rules allow certain collision coef- 
ficients, most notably Voooo> to become small in compar- 
ison with Vo2ii- In our numerical example we examine 
a situation where only the elastic collision coefficients 
Vbooo = Vo2ii/15 and Vbioi = V&211 are taken into ac- 
count. A further interesting feature of this system is the 
fact that the dominant collision mechanism and the atom 
cavity are realized with laser fields whose intensities can 
be varied independently from each other. This means 
that while the ratio Vijki /Vq2h is virtually fixed once the 
experimental setup has been optimized there is still con- 
siderable latitude in the choice of the relative magnitude 
of V0211 compared to Kq. 

For the numerical simulations involving the three- 
mode model to remain tractable, we chose parameters 
such that the laser linewidth is of the order of the natu- 
ral linewidth kq and the equilibrium laser mode popula- 
tion of the order of n = 50. With the help of Eq. ( g^ ) 
one estimates from the condition 033(2/^0) = 1 that 
Vqooo = 0.035k . This also fixes Vqioi = V0211 = 15Voooo- 
From Eq. ( p5| ) one obtains T = 5.63 x 10~ 3 ko as a suit- 
able value for the pumping parameter by again using the 
condition (733(2 / 'kq) = 1. Choosing k\ = 20kq yields 
N = 18.3. In this way all parameters which enter the 
two-mode model are determined. The relative orders of 
magnitude appear to be reasonable. The numerical cal- 
culation yields a half-life time of about 0.5/«o and no = 68 
which agrees reasonable well with our stated objective. 

In order to put this result into the right perspective 
it is necessary to perform corresponding calculations for 
the three-mode model also. This is because the two- 
mode scheme assumes that V0211, 1^2 - * 00 with the ratio 
r = 4V" 2 211 /k2 held fixed, without taking into account 
the need to assign finite values to the ratios Vijki /Vo2ii- 

One of the most important results of a systematic com- 
parison between the predictions of the two-mode and 
the three-mode models is that in the presence of elas- 
tic collisions the equilibrium populations may react as 
least as sensitively to the effects of "non-adiabaticity" 
as the linewidth. In other words, as the values of V0211 
and K2 are lowered with T fixed, uq can drop signifi- 
cantly below its adiabatic value even before the linewidth 
is considerably affected. This occurs as soon as the 
value of V0211 becomes of the order of the other colli- 
sion coefficients. This may appear surprising in view 
of the fact that in the two-mode system the popula- 
tion dynamics are completely independent of the elas- 
tic collisions. A possible explanation is that the levels 
|no, ni, 712) and |tiq + 1, n\ — 2, + 1) which are coupled 
by the coherent pumping experience an effective detun- 
ing due to the collisions. Hence, the predictions of the 



two-mode system for a given set of parameters should 
always be compared the full three-mode simulations to 
assess their reliability. In the present example, how- 
ever, the population decrease in the three-mode system 
(with K2 — 4Vo 2 2 n/r = 201ko) is not very pronounced 
(no = 48) while the linewidth hardly changes compared 
to its two-mode value. These results indicate that in the 
atom laser scheme of Refs. H|ll],[l2| it should be possible 
to obtain a linewidth of the order of or smaller than the 
natural linewidth k of the laser mode. 



IV. RELATIVE PHASE BETWEEN TWO ATOM 
LASERS 

The question of how a definite relative phase can be es- 
tablished between two independently created, interfering 
Bose condensates has been the subject of intensive the- 
oretical study in the recent past p6|-pl[. It was shown 
that one needs not resort to the concept of spontaneously 
broken symmetry to model the creation of interference 
patterns. Rather, these structures are brought about by 
the very act of observation: A sequence of measurements 
leads to the creation of an entangled state of the conden- 
sates pair which can be assigned a relative phase. How- 
ever, this phase varies randomly from one realization of 
the experiment to the other. 

In this section we examine whether such a relative 
phase is also created when the output beams of two in- 
dependent atom lasers are brought to interference. The 
purpose of this study is to determine to which extent the 
concept of measurements-induced phase built-up retains 
its validity also in open systems. The entanglement be- 
tween the atoms in the two laser modes is now influenced 
by the continuing pumping processes, and it is natural to 
ask if a well-defined relative phase can persist and what 
its properties are, e.g. with respect to diffusion. The 
question is also of interest with regard to related studies 
of optical systems |36| . 

We consider a gedanken experiment where the outputs 
of two two-mode atom lasers are brought to interference 
with the help of a 50:50 atomic beam splitter^. In the 
spirit of Ref. |50|, the effective non-hermitian Hamilto- 
nian describing the Schrddinger-like evolution of the two- 
laser system is 

H eff =H a + H b (26) 

where 



4 The use of a two-mode rather than a three-mode laser 
model simplifies the numerics while retaining the essential 
physical characteristics of the binary collision model. Fur- 
thermore, in the Monte Carlo wave function of a single two- 
mode atom laser a single number state is populated, so that 
it closely resembles the state of an isolated Bose condensate. 
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[K aja + «i(iV + l)a{ai 



t i r t 2 t 2i 



(27) 



is the effective Hamiltonian for the first atom laser. The 
Hamiltonian for the second laser is defined similarly, but 
using annihilation operators hi. The definition of the var- 
ious quantum jump operators is obvious from Eqs. j26| ) 
and ( p7|) except for those which describe atoms leaving 
the laser modes after passing through the beam splitter. 
They are given by 



(28) 



It is a simple matter to perform quantum Monte Carlo 
simulations within this framework. The general form of 
the simulation wave function can be written as 



!*(*)> = 



noi,m.i 

E 

™oi=noi,n 



Cn 01 (t) | noi,rc.H,no2 = N L - n i,ni 2 ), 



(29) 

where denotes the number of atoms in mode i of 
laser j and Nl is the total number of atoms in the two 
laser modes. The distribution function of the relative 
phase between the laser modes is conveniently expressed 
in terms of the overcomplete set of phase states [pCj 



i 



(aje i$ / 2 +ble-^/ 2 ) NL \0,0,0,0). 



(30) 



Because of the relation 

(oo + 6o)l*>JVx. - v / 2ivZcos($/2)|$) JVL (31) 

one can think of the state | $) m l as possessing a definite 
relative phase $ between the laser modes within the con- 
text of the interference experiment. The wave function 
can be expanded on these phase states as 



!*(*)) 



n d<S> 



c($,t) |$)jv £ 



(32) 



where the phase distribution amplitude c($,i) is given 
by 



c(*,i) = 2^/ 2 



"01, ma 

E 

"01=«01/ 



n Q1 l(N L - n 01 )\\ 1/2 



(33) 



Several physical processes influence the phase distri- 
bution: they are (i) the loss processes induced by the 
action of the operators c± leading to the entanglement 
of the atom lasers; (ii) the pumping into the laser modes 
described by the action of aj and 6 ; (iii) the collision 
interactions taking place between quantum jumps and 



described by H c ; and (iv) the non-Hermitian time evo- 
lution described by the T-term in Eq. (^?j). Whereas (i) 
and (iii) are also present in the case of Bose condensates, 
the influence of (ii) and (iv) occurs only for the laser. 
It is not a priori evident that the distribution function 
c(<I>, t) still yields a well-defined phase in the presence of 
these terms. 

Figure |] shows typical phase distributions |c(<I>)| 2 ob- 
tained from Monte Carlo wave function simulations of 
the two-laser problem. All examples in this section are 
calculated for the parameters T = 21.9ko, «i = 54.6ko 
and N = 2.2. Figure ||(a) illustrates the collisionless case 
Voooo = 0. Obviously, in this case the state of the system 
can be ascribed a well-defined phase. The distribution is 
symmetric with respect to $ = 0, reflecting the fact that 
the method of detection does not discriminate between 
the states |4>) and | — $). However, it is much broader 
than for the corresponding phase state \&)n l shown in 
Fig. |^(b) for comparison. That means that in general an 
approximation of the instantaneous Monte Carlo wave 
function by the superposition of two phase states is not 
appropriate. 

The influence of elastic collisions on the phase distri- 
bution is demonstrated in Fig. |](c). There, the value 
"Kjooo = 0.25ko was chosen, yielding a half-life time of 
T i/2 = 0.33/ko for the correlation function of the sin- 
gle atom laser, a value much smaller than the value 
T1/2 = 142/ko for the collisionless case. Figure ||(c) shows 
that the collisions tend to broaden the phase distribu- 
tion. Furthermore the symmetry of the distribution is 
lost with respect to both the height and the center of the 
two maxima. Nonetheless it is still reasonable to asso- 
ciate a well-defined phase with this state. 

From a formal point of view it is not too surprising 
that the Monte Carlo wave functions retain their definite- 
phase character in the presence of the pump processes. In 
the entangled state ( |29|) which is created by the interfer- 
ence at the beam splitter the coefficients c„ 01 (i) of signif- 
icant weight are concentrated in a small interval around 
noi = Nl/2. That means that the processes (ii) and (iv) 
change the wave function slightly, but without altering 
its basic character. 

Having established that a well-defined phase still exists 
in the presence of pump and collisions, we now turn to its 
dynamical properties, in particular its diffusive behavior. 
To this end we associate an instantaneous phase with the 
Monte Carlo wave function |\&(t)) via the relation 



cos 2 [$(i)/2] 



(*(t)\(at, + bl)((ao + bo)m)) 
2<¥(t)|aSao + 6£&o|¥(i)) 



(34) 



compare with Eq. (J3l|). Note that Eq. ( p4| ) accounts for 
the fact that the phase is only well-defined for the interval 
between and w. Note also that the use of other equiva- 
lent definitions of the relative phase, e. g. the maximum 
of the phase distribution |c($)| 2 , is possible. Eq. ( |3^ ) 
allows one to numerically evaluate averages of the form 



Hit) 



! [$(i)/2]) 



4>(t=0)=*- 



(35) 
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which measures the diffusion of an ensemble of Monte 
Carlo wave functions which all have the same initial rel- 
ative phase One can also determine the variance 



V^t) = (cos 4 [$(i)/2]) $(t=0)=i - A|(t), 



(36) 



where cos 4 [<£>(£)/2] is calculated as the square of 
cos 2 [$(£)/2] as inferred from Eq. A third quantity 

of interest is 

M(t) = (({cos 2 [$(i)/2] - cos 2 [$(i = 0)/2]} 2 )), (37) 

where the inner set of brackets indicates an average over 
quantum trajectories with the same initial $(£ = 0) and 
the outer one averaging over the different values of 4>(t = 
0). 

The quantities introduced in Eqs. ([35|) - (^) allow one 
to study the diffusion properties of the relative phase as 
defined by Eq. (|34|). If this definition is meaningful then 
there should be a connection between this result and the 
phase dynamics of an individual atom laser which, as we 
have seen in Sec. Ill, is well described within the frame- 
work of the linearized fluctuation approach. In partic- 
ular, the probability distribution for the phase 4>o of a 
single laser is given to a good approximation by 



i 



exp{-(</> o -0 o ) 2 /[2a 33 (i)]} (38) 



where </>o is the laser phase at t = and cr 33 (£) is given 
byEq.ll). 

Eq. (p8|) predicts that the time dependence of the quan- 
tities defined in Eqs. ( |35j ) - (|37j) are of the form 

Aj(t) = !{l + exp[-a 3 3(*)]cos$}, 

V*(t) = ±{l-2exp[-2a 33 (i)]cos 2 <l + 

exp[-4o- 33 (t)]cos2$}, 
M(t) = i{l-exp[-a 33 (t)]}. (39) 

where the calculation of M (t) assumes that all initial 
phase differences have equal probability. 

Fig. ^| compares these predictions to a Monte Carlo 
evaluation of Eqs. (35) - (p7|). Fig. ^|(a) and (b) show 



Ao(t), V (t) and V^/ 2 (t) for the collisionless case and 
the other parameters as in Fig. 5. Figure ||(c) depicts 
M(t) for Vbooo = 0.25ko- In all cases the behavior of 
the numerical results is well described by the expressions 
(|39|). A similar degree of agreement was also found in 
further examples. These results indicate that it is indeed 
possible to associate a definite relative phase with single 
simulations of the interference experiment and that this 
heuristically introduced phase behaves as predicted by 
the description of the atom laser in terms of stochastic 
processes Q. 

Experimentally, the phase can be observed (as soon as 
the correlation time is much larger than the mean time 
delay between atom emissions) by monitoring the ratio 



of the output intensities at the two beam splitter ports. 
Thus, the notion of spontaneously broken symmetry ap- 
pears not be the most natural concept to explain these 
experiments, either. 

As to the numerical simulations it should be remarked 
that <E>(i) as defined through Eq. (|34|) never assumes the 
values 0° or 180°. Instead, it always remains in the region 
between approximately 10° and 170°, for the collisionless 
case (and the specific parameters of our simulations), and 
a somewhat narrower interval in the presence of colli- 
sions. In fact, the numerical calculations of A^it) and 
Vo(t) were performed using the actual minimum value 
of <I> as initial value. For this value cos 2 ($/2) is still 
close to 1. At first sight, this observation might invali- 
date the description of the phase dynamics in terms of 
an unbounded stochastic process. However, this seems 
rather to be a shortcoming of the simple definition (|34|). 
It is more appropriate to think of the phase as having 
at every instant a well-localized but not infinitely nar- 
row distribution. Equation (|Z]) assigns a definite value 
to this distribution but may become inaccurate at the 
extreme values of cos 2 [$(t)/2]. In the calculations, this 
is compensated for by the fact that the probability that 
cos 2 [$(t)/2] assumes its actual minimum or maximum 
is somewhat larger than the probability for intermediate 
values. Furthermore, it should be noted that in Fig. ^|(c) 
the numerical result assumes quite exactly the predicted 
long-time value of 0.25. This justifies thinking of the 
probability distribution of the relative phase as being ap- 
proximately constant. 



V. SUMMARY AND CONCLUSIONS 



Using a linearized fluctuation analysis we have ob- 
tained analytical estimates for the influence of elastic col- 
lisions on the linewidth of a two-mode atom laser system. 
In case elastic interactions between laser mode atoms pro- 
vide the dominant collision mechanism the power spec- 
trum is Gaussian and the linewidth scales as Voooov^o 
for weak enough collisions. It undergoes a transition 
to a Lorentzian, with a linewidth proportional to no, 
as Voooo increases. A Schawlow-Townes-like behavior 
of the linewidth is recovered if collisions between pump 
and laser mode atoms are dominant. In that case the 
lincshape is Lorentzian with a linewidth proportional to 
Vqiqi/uq. The accuracy of the analytical predictions was 
confirmed by Monte Carlo wave function simulations. In 
addition, it was shown in a numerical study that the 
recently proposed atom laser scheme of Refs. 12 
might be capable of producing linewidths of the order 
of or smaller than the natural linewidth kq of the laser 
mode. In the course of the numerical study the useful- 
ness of the analytical estimates in the search for suitable 
operating regimes became apparent. 

In a second aspect of the phase dynamics, we inves- 
tigated the relative phase between two interfering atom 
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lasers. This study extends earlier work on the relative 
phase of Bose condensates to the realm of open systems. 
The analysis of suitably defined phase distribution func- 
tions showed that a physically meaningful relative phase 
can be ascribed to single runs of an interference experi- 
ment. A close connection with the phase dynamics of the 
single atom laser, as described by a Fokker-Planck equa- 
tion, was established by studying the diffusion properties 
of this relative phase. 
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FIG. 1. Schematic three-mode atom laser scheme. 



FIG. 2. Comparison between numerical calculation of the 
correlation function (full curve) and the result of Eq. jl^ ) 
(dashed) for parameter values P/ko = 15, ki/ko = 100, 
N = 1.85, and Vqooo/ko = 2.5, Vbioi = 0. 



FIG. 3. Inverse half-life time T x h of the correlation func- 
tion as a function of Voooo for parameter values P/kq = 15, 
«i/ko — 10, N = 20.3 and Voioi — 0. Shown are the results 
of Eq. @ (full curve), Eq. @ (dashed), Eq. Q (dotted), 
and numerical calculations (•). 







FIG. 4. Inverse half-life time r x ,\ of the correlation func- 
tion as a function of Vbioi f° r parameter values T/kq = 15, 
Ki/ko = 10, N = 20.3 and Vbooo = 0. Shown are the results 
of Eq. ( |l7| ) (full curve), Eq. ( pB| ) (dashed), and numerical cal- 
culations (•). 

FIG. 5. Examples of phase distributions |c(<S>) | 2 : (a) typical 
atom laser state for the collisionless case (other parameters as 
given in the text); (b) phase state as defined by Eq. (^o|) with 
$ = 45° and Nl = 120; (c) typical atom laser state in the 
presence of collisions (Vbooo = 0.25ko). 

FIG. 6. Comparison between numerical calculation (full 
curves, averages over approx. 250 simulations each) and an- 
alytical prediction (dashed) according to Eqs. (^) for: (a) 
Ao(t) with Voooo = 0; (b) V (t) and V n/2 (t) with Vbooo = 0; 
(c) M(t) with Vbooo = 0.25ko. Other parameters as given in 
the text. 
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